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by Peter M. Sockol 
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ABSTRACT 

An integral boundary layer procedure has been developed for the com- 
putation of: viscous and secondary flows along the annulus walls of an 
axial compressor. The procedure is an outgrowth and extension of the 
pitch-averaged methods of Mellor and Horlock. In the present work sec- 
ondary flow theory is used to develop approximations for the velocity 
profiles Inside a rotating blade row and for the blade force deficit 
terms in the momentum integral equations. The computer code based on 
this procedure has been iteratively coupled to a quasi-one-dlmensional 
model for the external inviscid flow. Computed results are compared with 
measurements in a compressor cascade. 


NOMENCLATURE 

profile parameters, eq. (9) 

i 

angular blade spacing, fig. 2 
axial projection of blade chord 
entrainment coefficient, eq. (19) 
skin friction coefficients, eq. (18) 
force deficits, eqs, (16) and (17) 
law of wall function, eq. (25) 
external and local values of blade force 
n power-law exponent, eq. (10) 

P,p external and local static pressure 

q ,q components of friction velocity 

r radius vector from axis 


A,B 

b 

C 


C c ,C,. 
fs ’ ft 

D. ,D 
t* x 

f(y c ) 

F .f 


Research Engineer, Computational Fluid Mechanics Branch. 
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s,t,y intrinsic coordinates, fig. 2 
U ,U„ components of wheel speed Q x r 

£3 L 

u rotational velocity field, w - w 

V,v external and local absolute velocities 

17, w external and local relative velocities 

W magnitude |w( of velocity K 

x,0,y cylindrical coordinates, fig. 1 
a arc tangent of end-wall slope, fig. 1 

8 external flow angle, fig. 2 

Ap pressure difference across pitch, 

8 shear layer thickness 

6' displacement thickness, eq. (14) 

- C {: power-law function, eq. (10) 

>1 y integral of vorticity component 

n* integral of n a , eq. (?) 

momentum thickness, eq. (14) 

K , K intrinsic curvatures, eq. (15) 

SC 

k ,k cylindrical curvatures 

x y 

X reciprocal length, ir/fcb cos 8) 

vorticity parameter, eq. (35) 

v thickness integral, eq. (23) 

ot 

p density 

T scaled 0 coordinate 

II Tq q wall shear-stress component 

<j» velocity potential of external flow, eq. (1) 

$ vector potential for velocity u, eq. (1) 




3 

S angular rotation vector 

u absolute vorticity 

Subscripts : 

0 denotes value at end wall, y « 0 

1,2 denotes value at blade surface 

Superscripts: 

denotes pitch-average 

+ 

Denotes average of blade surface values 


1. INTRODUCTION /• 

t! i; 

The three-dimensional viscous flows along the annulus walls of an 
axial compressor are known to play a major role in the performance of the 
machine. The turning of the shear layers by the blades and the relative 
motion between the end wall and the unshrouded blade tips produce secon- 
dary flows most of whose energy is ultimately lost. The displacement 
effects of the viscous layers can produce substantial local modifications 
in the pressure field of the external inviscid flow. The prediction of 
the significant features of these flows is the goal of the present work. 

Most of* the existing analyses of end-wall flows have followed one 
of two approaches. The first employs secondary flow theory in the small- 
shear large-disturbance approximation. An upstream primary vorticity is 
convected and distorted by a turning irrotatlonal flow. The resulting 
streamwise vorticity component Induces secondary velocities which are 
calculated by a stream function solution in a cross-sectional plane. 

Work in this area has been reviewed by Horlock and Lakshminarayana 
(ref. 1) and more recently by Horlock (ref. 2). This technique can give 
quite acceptable values for the distribution of flow angle downstream of 
a blade row, but it cannot predict the blockage and loss associated with 
the end-wall flow. 

The second approach, developed by Mellor and Wood (ref. 3) and 
Horlock and Perkins (ref. 4), employs integral boundary layer techniques. 
The three-dimensional equations of motion are averaged across the blade 
pitch and then integrated across the shear layer normal to the end wall. 
While the usual boundary-layer approximation is applied to the shear 
stress, the pressure variation normal to the end wall is not neglected 
since the blade spacing may be comparable to the shear layer thickness. 

In principle these equations can predict both the blockage and loss in 
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the end-wall region; however, in common with all integral procedures, 
many approximations are required to evaluate Important terms in the 
equations. One premising approach, which has been investigated by 
Horlock (refs. 4 and 5), is to use secondary flow theory to model some 
of the terms arising from three-dimensional effects. To date the in- 
tegral procedures have supplied useful predictions of the overall effects 
of end-wall flow o a blade-row performance. They have not, however, pro- 
vided the information on conditions within a blade passage needed in our 
expanded efforts to calculate the flow through a compressor stage. 

In the present work secondary flow theory is employed in a number 
of ways to improve the integral boundary layer technique. First the 
velocity profiles inside a rotating or stationary blade passage are ex- 
pressed in terms of Integrals of the pitch-averaged absolute vorticity 
components, which are in turn approximated by simple parametized forms. 

The velocity profiles are used not only to evaluate the various Integral 
thicknesses , : but also to calculate the effects of the pressure gradient 
normal to the end wall. Finally, the small-shear large-disturbance 
approximation is used to obtain the streamwise vorticity near the outer 
edge of the shear layer. Additional auxiliary information is provided 
by three-dimensional extensions to the law of the wall and to Head's 
entrainment equation (ref. 6) . 

The solution of the resulting ordinary differential equations for a 
multi-stage axial compressor has been developed into a computer code. 
Initial calculations with this code were performed for the case of the 
compressor cascade of Papalllou, Flot, and Matheu (ref. 7), The external 
flow was obtained from an inviscid calculation (ref. 8) in which the 
blockage was distributed approximately in accord with the measured values. 
The resulting end-wall predictions were extremely sensitive to the assumed 
blockage variation inside the blade passage; an Interactive calculation 
was apparently necessary. For simplicity in checking out the end-wall 
code a quasl-one-dimensional model was devised to represent the external 
inviscid flow (ref. 9). The axial velocity depends on the combined end- 
wall and blade blockage and the tangential velocity varies at a rate de- 
termined primarily by the turning of the blade mean camber line. Mien 
this model is allowed to interact iteratively with the end-wall calcula- 
tion, a satisfactory prediction scheme results. Comparisons between 
theory and experiment are / presented for the data of reference 7. 


2. BASIC ASSUMPTIONS 

The fluid equations are written for a reference frame fixed with 
respect to a rotating blade row. The two coordinate systems used in this 
frame are shown in figures 1 and 2. One is a cylindrical system with x 
along a surface meridian, 0 circumferential, and y normal to the end- 
wall. The second is an Intrinsic system with s and t, respectively, 
along and transverse to the streamlines of^ the external Inviscid flow. 

The variations of the external velocity W with y and flow angle 3 
with 0 and y are neglected. 
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In the evaluation of the secondary flow, and subsequently, the 
blade force deficit, the blade spacing rb Is taken to |ie of the same 
order as the shear layer thickness 6. Hence, we assume '.hat velocity 
gradients In the y and t directions are both large compared to those 
in the s direction with the possible exception of a small region near 
the blade leading edge. Both <$ and rb are assumed small compared to 
the blade height. 

In the vicinity of the blades we neglect the blade boundary layers 
and assume slip flow along each surface. Downstream of a blade row we 
neglect the viscous wukes and assume the secondary flow to be confined 
by vortex sheets which follow the external flow streamlines from the 
trailing edge. At the leading edge of a following blade row the pitch- 
averaged vorticity and velocity components are taken to be continuous. 
Although relative motion is permitted between the blades and the end- 
wall, no allowance is made for finite tip-clearance flows. 

Finally, we note that compressibility is included by allowing the 
external density p to vary with x. The variations of p with y 
and 0 are neglected. 


3. VELOCITY PROFILES AND SECONDARY FLOW 

In this section an approximate solution for secondary flow is used 
to develop expressions for the velocity profiles. The analysis is an 
extension of a similar one by Horlock (ref. 5) . A more detailed presen- 
tation is given in reference 9 . 

ji 

y ^The^veloclty^ w relative to) the blades is split into two parts, 
w - W — u , with W and u expressed as 

-v •+ 

W “ V$ - n x r , 

and u **- 0 for y » 6 with V • f = 0. 
the vector potential t are related by 

V ij; » ti) (2) 

We neglect Bernouli surface rotation and set M => tfj r » 0. Then the com- 
ponents of u can be written ^ ” 

3*. 9 

U s = " 3y~’ u t * 37”’ u y ' 3x~" ■ 7b sec B 

where r - 2(0 - Q m ) /b and the partial 3/3x is taken at constant t. 
Now 3/3x = cos 6 3/3s and we assume 3/3s = 0, hence equation (2) is 
approximated by 


u = ?x{ (1) 

The absolute vorticity m and 
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'a 2 1 4 ,2 r_' 

^y 2 V at 2 , 


£ S Id 


with X =» ir/(rb cos 8) . The approximate solution for i|* s is 

S I ^B ^ C0S I T 
* f ^ 

? B # n B <$e~ Xy - j n 8 (y x )dy x 




n s (y)dy 


(4) 


(5) 

( 6 ) 


(7) 


where t|) a is the average of tj/ across the pitch and 3n.,/3y 5 w a with 
n a = 0 for y > <5. The primary component tt I s simply approximated by 
its average value 




\ (y l )dy l 


( 8 ) 


Equations (2) and (5) to (8) can be combined to give velocity profiles 
once approximations for the vorticity components are established. 


Absolute Vorticity 


Consider first an axisyminetric boundary layer on a stationary end- 
wall upstream of a rotor. For this case X = 0 . The components n t t 
and ti 8 i in the absolute frame are approximated by 


with 


n t » “ -(V - v flf ) - -V(l - C) 


n 8 , “ ~v tf * -V(1 - c)cfA(i - O + BC] 



y — 5 
y > 6 


(9) 


CIO) 
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The parameters S, n. A, and- D will be permitted to vary with x. The 
streamwise component v^i has been given a simple power-law profile and 
the cross-flow v" t » approximates Johnston's triangular form /ref, 10) 
whenever A and 13 have the same sign. 


The components r) t and n a In the relative frame are given by 

n c » ~(w + u a )(i - ?) - u t (i - OcCACi - z) + m 
n a * -(w + u g )(i - OcEACl - c) + Brj + u t (i - O 


(li) 


where U a are components of the wheel speed fi * r. Inside the rotor 
X > 0 and W varies with t as well as x. We retain equations (11) 
but recognize that the velocity profiles obtained from them will not in 
genera.«- satisfy no slip at the wall. This is especially true at a rotat- 
ing rotor hub. The no slip condition is met by introducing a very thin 
region adjacent to y ** 0 in which the velocity profiles are matched to 
a three-dimensional law of the wall. The mechanics of this are taken up 
in the section on wall shear stress * 


Leading-Edge Region 

Whenever there is significant cross-flow relative to the blades up- 
stream of the leading edge, the assumption 3/Ss £ 0 Is invalid. The 
cross-flow must be turned by the blade, and it is only after some dis- 
tance inside the blade row that the passage vortex represented by equa- 
tion (6) is formed. This _situation can be modeled by the addition of a 
tern in x and y to y s in equation (6). For the case of a rotor 
tip, as described in the preceding section, this term has the value 
Hgd[l - exp(-Xy)] at the leading edge and it can be represented as a ! 
Fourier integral inside the blade row (ref, 9) . The displacement effect 
arising from this term can be shown to decay exponentially with distance 
x from the leading edget This might be used to obtain a suitable approx- 
imation to the Fourier integral. This possibility is being investigated. 


4. VISCOUS EQUATIONS 

The momentum Integral equations in the intrinsic system together 
with the entrainment equation are written In a form similar to that of 
Hlrsch (ref. 11), A complete derivation is given in reference 9, In 
this section the preceding velocity profiles are used to develop approxi- 
mations to the blade force deficit, wall shear stress, and entrainment 
rate. 

The momentum integral equations are 


(19 


df + % £ 1 "< cb(,w2) + 


i / .* aw , * 


aw 

at 


- ie.9. + K 0 fcfc + 2k U6. s D , cots 

t ts S tt X t X 




( 12 ) 


d0 


¥ + 8 tX £ b P« 2 > %?«« tS + *t(> 


2k x 5 < * -»x ai “ + D t + k C f t 


(13) 


where 


W 6 
a a 




w 2 e 


<*8 



(W ~ w )dy 
a a ' 


<»„ - Vs d » “ ,w X ■ 


u dy 
a 7 





u u a dy 
a S 


/ (U) 


as . 

k = - tt- + k cos B , 
s at |, x * 


K t " aiT - K X sin B 


(15) 


D » -(pW z ) 


x 


{ S '. -1 JL 

dx 



(P - p)dy o -(pw2) 


-1 d 



3P 


d x / 3y 

J 0 


y dy (16) 


D t = (pW 2 ) 


-1 



(F t - f fc )dy « (PWV 1 



3f t 

ay" 


y dy (17) 


°Jo ■ 2 ^ 52 )' 1 v 


(18) 


with k = -3r/r 3x, U = U/W, and f » -AP/(rb cos 6). Here AP - P - P 

X C it X 

is the pressure difference across the pitch. The upper limit on y is 
taken very large because the secondary flow extends well beyond the edge 
of the shear layer, y « 6. 


The entrainment equation is written 
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(6 cos g - <S* + fi“) + (6 cob g - 5* + <$“) ~ ln(rbpw) “ Cg (19) 
where C„W is the entrainment rate at y » 6 and 



The integral thieknesbes are evaluated from the profiles of equa- 
tions (3), (6), (8), and (11) in reference 9. Th e ave rages in equa- 
tions (12) and (13) are approximated by setting w a wg * w a wg, Estimates 
of the neglected terms indicate that they contribute less than ten per- 
cent even for highly loaded blades. 


Force Deficits 

Let 3/3s be of order one, 0(1), and assume that rb is 0(5) , 
then from the y momentum equation 

™ 3 £ 

(PK 2 )' 1 f- ■ 0(1) . (pi? 2 )" 1 - o(i) 


When these relations are substituted into equations (16) and (17), it is 
seen that D is 0(5 2 ) and Dj. is 0(6). Since equations (12) and (13) 
are only correct to 0(6) we neglect D x . An equation for fj. is ob- 
tained by evaluating t! a y momentum equation at the blade surfaces 
where w^ «=> Q. 


5 

3y 


w 

s 

rb 


JL 

3x 


(Aw y ) 


rb" 


+ 

+ A 

3x rb 


. 1 . 

by 


Cw,. 


Aw y )sec 


+ 2tt (wj + U) 



tan 0 + viscous 


terms 

( 21 ) 


4 * 1 

where W Q = j (w al 4- w^,,) and K y = -3r/r 3y. The integrand in equa- 
tion (1/) is only significant for y of 0(6) or larger where 
Wg 1 * W + = W and the viscous terms in equation (21) are small. Also from 
equations (3) and (5), away from the leading edge, w y * 0. The first 
term on the right of equation (21) must be retained /however , since 
AWy = Wy2 - Wy^ is 0(w t ). With these approximations, substitution of 
equation (21) Into equation (17) gives 





( 22 ). 
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dV v j _ 

D t s„~co8 0 ^ M cos 0 ^ ln(rbW cos .A) 


- 2Ky(bin 0 + U)v g sin 


with 


Hv * (rb cos 0) 

U 


- 1 / 


ACW a - w a )y dy - ~ 



Au a y dy (23) 


The thicknesses v are evaluated in reference 9. 
a 


Wall Shear Stress 

We assume the flow in a thin region adjacent to y “ 0 is collateral 
and obeys, a three-dimensional law ?>£ the wall, 

* s " w s0 3 7 t “ w t0 " V <y +> < 24) 


where 


f (y + ) = o^r ln y+ + 5 > y + - yq/v (25) 

and w a Q are the components of the end-wall velocity relative to the 
blades. The skin friction coefficients ajfe related to the friction ve- 
locities q a by C£„ = 2qq a /w2. At y 0 , equations (3), (6), (8), and 
(11) give 


W B - “so ■ - ( “ U + V .0> + [<« + V - V 1 * 
w t “ w to ” - ^ u t + w to> + CW)nJ + [(w + u a )A + u t ]£ 

Matching equations (24) and (26) at y / & - y gives 


(26) 


q s £(y + ) • -<D S + w s0 ) + KW + U s ) - U t A]y 1/n 


q t f(y + ) * -(U t + w tQ ) + (X6)n* + t(w + u s )a + u t ]£ 1/n ^ 


t(r 


(27) 


with y, » y6q/v. Equations (27) can be solved iteratively for q, and 
q t . This formulation enables the flow to respond to sudden changes in 
the end-wall velocity such as those which occur at a rotor hub. Note 
that y 53 0.1 gives a good value for Cf in a two-dimensional boundary 
layer. 
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Entrainment* Rat- 


In two-dimensional flows the entrainment rate CgW Is given by 
(ref. 12) 



1 

O 


/y=6 


(^x\ 

Vw Vs 


C28) 


where \>‘j* is the eddy viscosity. Since Cg is n scalar, we assume with 
Brudshaw (ref. 13) that the same relation holds in three-dimensional flow. 
If we approximate v.j. by 


y * 

• *!» 



(29) 


where % is the mixing length, then equations (11) and (29) as y *> 6 
give 


v. r « (v, r ) iD V(l + B 2 ) 1/2 , 


(v t ) 2d * 


,2 JV 

M ny 


(30) 


with V ° V/W. The expression for (vq^Q is not accuruce enough to be 

used in equution (28). Hence, its derivative is replaced by Head’s empir- 
ical relation (ref. 6). Thus, we approximate Cg by 


with 


C E * 0.03G6(H* - 3)~°’ 653 vU + B 2 ) 1/2 
11* “ 1.535(11 - 0.7)' 2,715 + 3.3 


(31) 

(32) 


H a 1 + 2/n 


(33) 


Equution (33) for the streamwise shape fuctor H follows from the power- 
law relation equation (10) , 


5. VORTICm TRANSPORT 


In the outer part of the shear layer, near y a 5, we assume that 
streamwise vorticity is convected by the external flow with only the 

primary component directly affected by viscous diffusion. This model 

is similar to that or Louis (ref. 14). In addition pressure and density 
are taken to be the same as in the external flow where P a P(p). For 
these conditions the transport equation for w s is (ref. 9) 


pW 


■9 

3s 


/M 

i \T/-* 


. w- 


3W 

3s 


2w t (« t W - ^ X U) 


(34) 
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Equation (34) it! evaluated with average values for bs a and W, From 
equations (11), as y -*■ 6 and g -* 1, 

s V^ ny C35) 

l» t a (1 + U H ) + U C B, » (1 + U S )B - U fc (36) 

with U a a U tt /W. Note that y must be much less than 6 before 5 de- 
parts significantly from one. Finally, equations (34) to (36) are com- 
bined to give 


7 $ < 37 > 

6. INTEGRATION OF THE GOVERNING EQUATIONS 

The integral and vorticity equations are written in a combined im- 
plicit fora 


37 h “ V 1 “ 1 e ° 4 «« 

The overburs are dropped and all quantities are understood to take their 
average values. From equations (12), (13), (19), (22), and (37) we write 


F^ “ S cos 


- § + 6 , 
X X 


f„ « e 
2 ax 


F 3 “ 0 tx + COB g » 


F 4 - p B /n 


( 39 ) 
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=1 ■ - F 1 & + C E 


°2 ■ - e sx H ln(cbpw2> - { x ff f + *l\s - “At 


- Vt + I C fa 


°3 “ Ax S' ln < rb ' iw2 > + "As - K t< 6 ss + 0 + *x US s 


> (40) 


- V y cos 8 ~~ ln(rbW) - 2K y (sin 2 B + U a )\> s + ~ C £t 


G 


a !® + 2 sec 3 “ (k, - k U) 

4 n p dx n't x 


P Numerical Procedures 

Equations (38) are rewritten in an explicit matrix form 


dY 


M(x,Y) ~ » G(x,Y) - F k (x,Y) 


(41) 


or 

s • M ' 1(G - V 2 D < 42) 

where Y is the unknown vector (<S,n,A,B), M is the Jacobian matrix 
9F/3Y., and F x = 3F/9x, The elements of M and F y are given in ref- 
erence 9 . Equation (42) is integrated by the fourtn-order Runge-Kutta 
method with Gaussian elimination used to find D at each step. 


7. DISCUSSION OF RESULTS 

The end-wall theory has been used to predict the low speed flow 
through a linear cascade of compressor blades as reported in reference 7. 
A detailed set of data for the heavily loaded case B has been obtained 
from the authors (ref. 15). In this flow a set of NACA 65-(12A10)-10 
blades at 30 degrees stagger was used to turn the flow from B = 56.6 
degrees to S = 25.2 degrees. The data taken upstream of the cascade re- 
veal the presence of a significant cross-flow extending to mid-span. The 


14 


S shear layer, however, la confined to a thickness of less than one fourth 
the blade height. Since the theory, In its present form, is incapable 
of properly treating this situation, the computation was performed with 
the upstream cross-flow set to zero. 

The comparison between theory and experiment is shown in figures 3 
to 6. Figures 3 and 4 show the tangent of the flow angle 3 and the 
axial velocity V x against dimensionless distance x/C a , where C a is 
the axial projection of the blade chord, as obtained from the quasi-one- 
dlmenslonal model of the lnviscid flow. The turning rate inside the 
blade row has been adjusted to produce the agreement in tan 3 shown 
in figure 3. However, the figure also shows that the real flow undergoes 
significant turning before it encounters the leading edge at x = 0. The 
under prediction of V x in figure 4 indicates that the real flow en- 
counters somewhat greater blockage than predicted in the computation. 
Figure 5 shows the streamwise displacement and momentum thicknesses , ^ 
6g/C a and 6 BB /C a , against x/C a . While the close agreement for S g 
inside the blade row is probably fortuitous, the general trends for both 
curves are highly encouraging. There is no point in Comparing cross-flow 
thicknesses, since the upstream cross-flow has been neglected in the com- 
putation. Finally, figure-. 6 shows the blade force deficit D t against 
x/C a . Although there is qualitative agreement, the predicted value is 
too small by more than a factor of two. The inclusion of a leading edge 
model which can properly account for the rapid turning of the upstream 
cross flow should go a long way toward removing this discrepancy. 


8. CONCLUSIONS AND RECOMMENDATIONS 

In its present form the end-wall prediction scheme has been shown to 
give at least qualitatively correct results for the pitch -averaged flow 
through a linear cascade of compressor blades. Naturally many more com- 
parisons with internal flow measurements are needed. However, at least 
three things must be added before the theory can be used to predict the 
end-wall flows in an axial compressor: 

1. A leading edge model which can treat the rapid turning of an up- 
stream cross-flow and its contribution to the blade force; 

2. A model for tip clearance flows and their effect on boundary layer 
growth ; 

3. Coupling of the end-wall prediction to a multi-stage though-flow 
calculation for the lnviscid flow. 
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Figure 1, - Ertd-all gsfimetry and cylindrical coordinate system. 



Figure 2. - Cascade geometry and intrinsic coordinate system. 




